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ABSTRACT 

Magnetic reconnection is believed to be the prime mechanism to trigger solar flares and accelerate electrons 
up to energies of MeV. In the classical two-dimensional reconnection model, the separation motion of chro- 
mospheric ribbons manifests the successive reconnection that takes place higher up in the corona. Meanwhile, 
downward traveling energetic electrons bombard the dense chromosphere and create hard X-ray (HXR) emis- 
sions, which provide a valuable diagnostic of electron acceleration. Analyses of ribbon dynamics and HXR 
spectrum have been carried out separately. In this Letter, we report a study of the comparison of reconnection 
electric field measured from ribbon motion and hardness (spectral index) of X-ray emission derived from X-ray 
spectrum. Our survey of the maximum average reconnection electric field and the minimum overall spectral 
index for 13 two-ribbon flares show that they are strongly anti-correlated. The former is also strongly correlated 
with flare magnitude measured using the peak flux of soft X-ray emissions. These provide strong support for 
electron acceleration models based on the electric field generated at reconnecting current sheet during flares. 
Subject headings: Sun: flares — Sun: magnetic fields — Sun: X-rays, gamma rays 



1. INTRODUCTION 

The "ribbon" structures of solar flares have long been ob- 
served at chromospheric lines (e.g. Ha at 656.3 nm and 
Ca II H at 396.8 nm). A ribbon in one magnetic polarity re- 
gion has its counterpart in the other magnetic polarity region, 
and both run parallel to the magnetic polarity inversion line 
lying between them. The well-observed separation motion 
of ribbons is a direct mapping of energy release via mag- 
netic r econnection in the cor ona reconnecting current sheet 
(RCS) (iPriest & Forbesll2002l) . the rate of which can be eval- 
uated as = (d/dt) J B da, where B is the corresponding 
magnetic field component perpendicular to the surface ribbon 
element da . Under a simplified two-dimensional reconnec- 
tion model (iForbes & PriestlfT984l) . the electric potential drop, 
y = J E dl, along coronal separator with length element dl 

equals the reconnection rate (p. By considering that the coro- 
nal separator has the same length as flare ribbon, the electric 
field can thus be expressed as £ = uB, where u is the ribbon ex- 
pansion velocity. This enables the determination of the elec- 
tric field at the RCS through observable quantities in the lower 
atmosphere, while large uncertainties may be induced in eval- 
uating u especially for ribbons with irregular shape. Alterna- 
tively, we infer the average electric field used in the present 
study as 
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where L is the length of ribbon, and the "+" and "-" sub- 
script denote the corresponding physical quantity measured 
for the ribbons located in the positive and negative magnetic 
fields, respectively. In this way, the reconnection rate </> can 
be calcul ated essentially based on progression of ribb on inten- 
sity (e.g.. lOiu & Yurchyshynll2005l:ISaba et"ani2006h . Impor- 
tantly, reconnecting magnetic fluxes from positive and nega- 
tive magnetic fields should be identical in principle; however, 
their actual measurements do n ot always yield a good balance 
(e.g.. iFletcher & Hudsonll200TI) . Therefore, the average value 
of E, as defined above, may be able to represent the overall 



strength of reconnection electric field at the RCS, which is re- 
garded as conceptually the most straightforward mechanism 
to dir ectly accelerate electrons to high energies (iLitvinenkol 
fT996h . 

A powerful diagnostic of accelerated energetic elec- 
trons produced by flares is the hard X-ray (HXR) emis- 
sion, the spectrum of which often appears as a power- 
law distribution in photon ene rgy [/(e) oc e"'''] (e.g., 
iTandberg-Hanssen & Emsli3ll988h . The hardness of X-ray 
emission, i.e., the power-law spectral index 7, implies a char- 
acteristic in the energy distribution of the electron flux bom- 
barding the dense target under the bremsstrahlung emission 
mechanism (lBrownlll971h . Specifically, a harder X-ray spec- 
trum (with smaller spectral index) indicates that more elec- 
trons are accelerated to higher energies. For this study, X- 
ray spectrum is measured for entire flaring region on the solar 
disc, thus the derived spectral index reflects the overall hard- 
ness of X-ray emission (thereafter referred as (7)) created by 
downward traveling accelerated electrons impinging on the 
chromosphere. 

Previous studies of temporal evolution of reconnection 
electric field E and HXR spectral index 7 in individual 
flares show that they are correlated and an ti-correlated, re- 
spectively, wifli that of the HXR flux (e.g.. lOiu et all 120041: 
iGrigis & Benzl2064b . which naturally suggests an anticorrela- 
tion relationship between E and 7. This is also hinted by par- 
ticle simulations with pre s cribed electric a nd magnetic fields 
(IWood & NeukirchI 120051: iLiu et alJl2009l) . which show that 
the electron energy spectrum hardens when the electric field 
increases. A direct comparison of E and 7 using observational 
data is thus needed to bridge the research of magnetic recon- 
nection and electron acceleration, which is believed to have 
an intrinsic causal relationship. 

In retrospect, only a single such attempt has been reported. 
It was the 2005 May 13 M 8.0 flare with ribbon-like HXR 
emissions (Liu et al.l |2007) . where spatially resolved E and 
7 determined by tracing each position along ribbon and imag- 
ing sp ectroscopy for e ach pixel, resp ectively, can be com- 
pared ( iJing et al.ll2007l:lLiu et al.ll2008l) . As ribbon-like HXR 
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Fig. 1 . — Calculation of reconnection electric field (E) from ribbon motion and derivation of hardness of X-ray emission (7) from spectral analysis for the 
2005 January 15 X2.6 flare. The general flare morphology is shown in (a), where an Ha image near flare peak time is overplotted with X-ray sources. The 
RHESSI images were integrated for 1 minute centered on the Ho image, and were reconstructed with the CLEAN algorithm using the front segments of detectors 
3-9 giving a FWHM resolution of ~9.8". The ribbon motion tracked using a central brightness cutoff value (see §|2) is depicted in (b), overplotted on a preflare 
SOHO MDI magnetogram. The X-ray spectrum (c) made in the same time interval as the X-ray images shows the model fitting using a themial Maxwellian and 
a nonthermal double power-law components, with background flux illustrated (grey). Temporal evolution of (E), (7), and RHESSI 25-50 keV photon flux are 
compared in (d). 



sources are not often observed, in this Letter, we take a dif- 
ferent approach and investigate the relationship between (E) 
and (7) using a sample of 13 two-ribbon flares. 

2. OBSERVATIONS AND MEASUREMENTS 

The chosen flares were all well observed on the disk 
with an orientation cosine factor of ^^0.7 (^46° from the 
disc center), exhibiting clear motion of ribbon separation 
with nearly full coverage of X-ray observations from the 
Reuven R amaty Hi^h En ergy Solar Spectroscopic Imager 
(RHESSI; iLin et alj|2002h . Ribbon observations at Ha and 
Ca n H spectral line were acquired by the Global Ha Network 
(ISteinegger et al. 2000) (10 events), the Solar Optical Tele- 
scope (SOT; ITsuneta et al . 2008) on Hinode (2 events), and 
the USAF/Optical Solar Patrol Netwo rk (OSPAN, formerly 
known as ISOON: iNeidig et alJfT998h (1 event). These im- 
ages have a cadence of 0.5-2 minutes and a pixel resolution 
of ^0. 1-1". For the sake of consistency, we use photospheric 
magnetic fields of the flaring region (^2" pix el"') measured 
with the Michelson Doppler Imager (MDI; IScherrer et alj 
Il995h on the Solar and Heliospheric Observatory {SOHO) for 



all the events. 

We take the 2005 January 15 X2.6 flare as an example 
to demonstrate our method of data reduction. Figure [TJ? 
shows an Ha image near the flare peak time superimposed 
with X-ray sources. This is a typical two-ribbon flare, with 
part of Ha ribbons occupied by HXR emitting sources (25- 
50 keV) and a lower energy X-ray source (12-25 keV) lying 
in-between. We developed a standardized procedure to derive 
the reconnection electric field from ribbon motion. Modules 
of this procedure include accurate image alignment, appropri- 
ate image destretching for ground-based observation, ribbon 
tracking using th e intensity-based binary masks method (e.g., 
'Saba et al1l2006i) . ribbon length estimation, and registering 
Ha images on the corresponding magnetogram by multiple 
feature matching. We consider uncertainties in this method 
that mainly stem from the following aspects: (1) The choice 
of the cutoff value that defines the brightened ribbons. We 
first determine a central cutoff value at the evolving edge of 
the ribbon based on the difference images relative to a fixed 
preflare frame. We then vary this central cutoff value within 
±20% and use 1 1 thresholds in total to characterize the bright- 



RECONNECTION ELECTRIC FIELD AND HARDNESS OF X-RAY EMISSION 



3 



TABLE 1 

INFORMATION OF FLARES IN THIS STUDY 



NOAA Location Cosine GOES GOES 1-8 A Maximum Electric Minimum HXR Time Difference' 
Event Date AR (degree) Factor Level Peak (UT) Field (£) (V cm"') Spectral Index (7) (minutes) 



(1)2002 Feb 18 


09830 


S20E16 


0.9 


Ml.O 


21:15 


2.3 ± 0.7 


6.84 





(2) 2002 Feb 20 


09830 


S18 Wll 


0.93 


M2.4 


21:07 


5.8 ± 1.3 


3.64 





(3) 2003 Apr 23'' 


10338 


N20 W22 


0.87 


M2.0 


15:56 


3.2 ±0.8 


5.1 


-1 


(4) 2003 May 27 


10365 


S07 W17 


0.95 


X1.3 


23:07 


7.7 ± 1.8 


2.7 


OF 


(5) 2003 Oct 29'' 


10486 


S15 W02 


0.96 


XIO.O 


20:49 


16.7 ± 5.9 


1.88 


1 


(6) 2004 Mar 30 


10581 


S05 EOl 


0.99 


C2.0 


23:08 


1.6 ±0.4 


6.94 





(7) 2004 Nov 4'' 


10696 


N08 E18 


0.94 


M5.4 


23:09 


6.2 ± 1.3 


2.98 





(8) 2005 Jan 15 


10720 


N15 W05 


0.96 


X2.6 


23:02 


9.2 ± 2.1 


2.46 





(9) 2005 May 13 


10759 


N12E12 


0.96 


M8.0 


16:57 


5.7 ± 1.4 


3.21 


-1.5 


(10) 2005 May 26'' 


10767 


S06E13 


0.97 


C8.6 


21:39 


2.1 ±0.5 


6.06 


1 


(11)2006 Jul 6 


10898 


Sll W32 


0.83 


M2.5 


08:36 


4.9 ± 1.7 


3.46 


-3 


(12) 2006 Dec 14 


10930 


S06 W46 


0.69 


X1.5 


22:15 


9.2 ±2.1 


4.0 


-1 


(13) 2007 Jul 10 


10963 


S07 E45 


0.7 


C5.2 


17:53 


2.1 ±0.5 


7.73 






Note. — Flare ribbon images of the events 1^ and 6-11, 5, and 12-13 were obtained with GHN, OSPAN. and SOT, respectively. 
" Time difference = time of maximum average electric field — time of minimum overall HXR spectral index.'' We only extract results from the 
later phase of these events, when we consider there exhibits classical separation motion of ribbons.*-' RHESSI most probably missed the main HXR 
peak of this event. 



ened ribbons at each time interval, in order to have a better 
evaluation of the reconnection rate; (2) Magnetic field mea- 
surement. MDI level L8 full-disk magnetograms are used, 
and a largest noise level of 30 G is taken into account. Since 
MDI magnetogram coul d be affected by h igh energy emis- 
sions during flares (e.g., lOiu & G^l2003h . we use magne- 
tograms just before the flares for calculation. We note that it 
is known that MDI data could suffer from Zeeman saturation, 
and the level L5 data coul d underestimate the flux density 
by a factor of ~0.64-0.69 (iBerger & Litesll2003h . However, 
in most of our studied events, the regions swept by flare rib- 
bons are not obviously affected by Zeeman saturation, and 
the values of level L8 calibrated field density are '-^L6 times 
those of the uncalibrated level L5 data'. Therefore, we do not 
expect the inherent limitations in MDI measurement to sig- 
nificantly alter out results; and (3) Ribbon length estimation. 
We measure the length of the brightened ribbons as defined 
above. While as the shape of flare ribbons may not be reg- 
ular and actual flare process may take place in a more com- 
plicated three-dimensional structure (see more discussion in 
§ nil, we take an error of 30% in ribbon length measurement. 
Moreover, caution has been taken when there is a complicated 
ribbon dynamics, only later stage of which exhibits classical 
separation motion and can be regarded as complying with the 
standard two-dimensional model (e.g., the 2003 October 29 
XIO flare; e.g.. lLiu et al. 2006). Thus in such cases (indicated 
in Table[TJ, we extract results only from the clear ribbon sepa- 
ration phase. Figure[TJ' depicts front of the X2.6 flare ribbons 
as color-coded lines using a central brightness cutoff value, 
which shows clearly the ribbon separation motion in opposite 
magnetic fields. The computed (E) (red) is plotted in Fig- 
ure [T}/. Its peak (~9.2 V cm"') is co-temporal with that of 
the HXR emission in 25-50 keV (black), which is pre viously 
found in other events (Oiu et al. 2004; Jing et al. 2005). 

For analysis of X-ray spectrum, we used the Object Spectral 
Executive (OSPEX) software package and applied the pileup 
correction, if necessary, in SolarSoftWare (SSW). We set the 
energy bin width to 1/3 keV, the standard finest binning, and 
choose the length of fitting time interval to be the same as 
the cadence of ribbon images centering on the time of each 

' |http : //soi ■ s t an f or d. edu /magnetic /Levi ■ 8/| 



frame, except in the case that we have to avoid the time when 
RHESSI attenuators change status. We model the thermal 
plasma and high-energy emissions with an isothermal com- 
ponent and nonthermal power-law distributions, respectively. 
Quite often, the nonthermal emission is best fitted using a 
double power-law function with two different spectral indices. 
In such a case, we take the spectral index reflecting the hard- 
ness of photons with medium energies (~20-50 keV) as the 
result for that event (see discussion in §|4]i. The uncertainty in 
the values of HXR spectral index may arise from the choice 
of the length of fitting interval and spectrum background, and 
is estimated to be r^5%. The X-ray spectral fitting results of 
the X2.6 flare around its peak time is shown in Figure [!}:, 
where (7) reaches ^2A6 in ~20-60 keV. The time profile of 
(7) (blue) is plotted in Figure[T]i, in comparison with those of 
{E) and the 25-50 keV photon flux. As usual, the minimum 
overall spectral index occurs when the HXR peaks, which is 
believed to be an intrinsic fea ture of electron accel eration in 
the flare impulsive phase (e.g.. lGrigis & Ben3l2004l) . 

3. RESULTS 

Following similar procedure, we calculate reconnection 
electric field and HXR spectral index for other events, and 
give the results in Table [T] The maximum {E) and the mini- 
mum (7) lie in the range of ^^1-20 V cm"' and ~2-8, respec- 
tively, which is generally comparable with w hat previously 
reported (e.g.. lJing et alll2005l; iLin et al.ll2001l) . In a majority 
of events, the times when {E) and (7) reach the maximum 
and minimum values, respectively, agree within one minute. 

Figure|2](to/:') shows the scatter plot of these two quantities, 
with a value of logarithmic correlation coefficient (CC) com- 
puted to be -0.94. That is to say, a larger electric field in the 
RCS corresponds to a harder X-ray spectrum (with smaller 
7) and presumably to a harder electron precipitation spec- 
trum, which is qualitatively consistent with the trend predicted 
by nu merical simulations ( Wood & Neukirch ,2005; Liu et all 
120091) . This high correlation relationship thus strongly sup- 
port the hypothesis that direct acceleration by the electric field 
generated by magnetic reconnection may play an important 
role in producing energetic electrons in flares. We note that 
X, M, and C-class flares appear to be three populations, which 
indicates that flares with larger magnitude are almost always 
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Fig. 2. — Scatter plots of the minimum overall HXR .spectral index (top) 
and the peak flux in 1-8 A (bottom) vs. the maximum average electric field 
in the RCS in a logarithmic scale. The dashed line shows the least-squares 
Hnear fit to the da ta points. The rms measurement error for the GOES 1-8 A 
channel is 16.2% )GarcidfT994h . 



associated with larger electric fields and harder HXR spectra. 

This tendency is more clearly seen in Figure [2] (bot- 
tom), where we make scatter plot of the maximum {E) ver- 
sus the peak flux in 1-8 A soft X-ray measured by the 
Geostationary Operational Environmental Satellite (GOES; 
[Kahler & Kreplin 1991). The logarithmic CC also has a high 
coefficient of 0.97 for all the data sets. This further im- 
plies that larger electric field strength may be responsible for 
stronger electron acceleration and hence stronger high-energy 
emissions seen in flares with larger magnitude. 

4. SUMMARY AND DISCUSSION 

Joining the research of magnetic reconnection and electron 
acceleration, we have determined the average electric field 
(E) in the RCS and the overall hardness of X-ray emission 
(7) from flare observations of 13 events. We then find a 
strong correlation relationship between both the minimum (7) 
and the flare magnitude corresponding to the maximum {E). 
This is, to our knowledge, the clearest direct observational 
evidence obtained so far of electron acceleration by recon- 
nection electric field in solar flares, whi ch has long been ex - 
plored but mostly by theoretical models (lAschwandenll2005h . 
We discuss the present results with related observation and 
simulation works. 

The present research is closest to that of iLiu et aT ( 2008"), 
who compared spatially resolved electric field and HXR spec- 
tral index for a single event. In their result, an anticorrela- 



tion relationship between E and 7 is also found, while for 
the E with similar strength, the corresponding 7 generally has 
smaller values than that in our study. A possible explanation 
for this discrepancy is that different from the average E and 
the overall 7 herein presented, they explored the local elec- 
tric field by tracing motion of every position along ribbon, 
and local spectral index in each pixel by applying an imaging 
spectroscopy technique. A larger error could be induced in 
calculating such E and 7, but it is more possible that differ- 
ent approach, i.e., global or local, matters in studying the flare 
energy release. In another study, jNind (I2008J found that 
is anti-correlated with 7 in two flares. However, one of their 
event is close to limb (with a cosine factor of ^^0.49), and they 
used ribbon observations in 1600 A that forms in upper chro- 
mosphere and transition region. More importantly, instead of 
(p, E is believed to be the physical quantity directly associated 
with electron acceleration. 

Previous studies of reconnection electric field usually use 
the form, E = uB, while u is not easy to be measured accu- 
rately (cf. iQuetaDEOOl. We use another form, E = d/L, 
where (j) ca n be evaluated in detail using ribbon intensity (e.g., 
ISaba et aT]l2006.) . Here we must note that the ribbon length L 
is an invariable in the standard two-dimensional model; how- 
ever, it usually involves rapidly during actual flares especially 
in the early phase. Although chromospheric ribbons may 
be regarded as an instantaneous response of energy release 
in the coronal current sheet through precipitation of acceler- 
ated electrons, recent observational studies have suggested the 
requirement of consid eration o f three-dirnensio nal magnetic 
structures (e.g.. lLiue t al. 20071: iJing et alj|2007h . The devia- 
tion from the standard two-dimensional model is also implied 
by the asymmetry of flare geometry. For the studied events, 
the ratio of reconnection electric field and ribbon length mea- 
sured in opposite magnetic field regions, i.e., ((p+Ll^)/((p-LZ^) 
and L+/L-, generally falls into the range of ^0.2-2 and 0.3- 
3, respectively. Without an applicable three-dimensional flare 
model, we thus use the equation in § [T] as an approximation, 
where we measure the length of ribbon when it fully devel- 
oped and treat it as a constant. This simplification may be 
justified, because most of our events do not show a significant 
change of ribbon length at times when the maximum {E) and 
the minimum (7) occur Using the time-dependent change of 
ribbon length as a probe to the magn etic reconnection p rocess 
is an area of on-going research (e.g.. lLee & Garvl[2008l) . 

Our observational results demonstrate that there exists a re- 
markable anticorrelation relationship bet ween E and 7. This 
qualitatively agr ees with model results (IWood & NeukirchI 
120051; iLiu et all 12009 ). in that both indicate a hardening 
of electron energy distribution with increasing electric field 
strength. Quantitatively comparison is beyond the scope of 
this study, as further assumptions need to be made before nu- 
merical results of electron spectra can be compared with ob- 
served p hoton spectra. A n important issue is that in the simu- 
lation of lLiu et al] (|2009|) . it is found that the accelerated elec- 
trons exhibit not only a power-law energy spectrum but also 
an exponential tail at the high energy end. The latter indicates 
that diffuse shock acceleration may be at work. Therefore 
we take only the spectral indices of the HXRs with medium 
energies as the results, when the spectrum can not be fitted 
using a single power-law function. This also implies that dif- 
ferent acceleration mechanisms may all contribute to electron 
acceleration during flares (Aschwanden 2005). Nevertheless, 
the presence of such clear correlations of observational signa- 
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tures provides strong support for electron acceleration models 
based on electric field generated at the RCS. 
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